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Abstract 

By using the Renormalization Group Equations in Chiral Perturbation The- 
ory, one can calculate the double chiral logs that appear at two loops in any 
matrix element. We calculate them in the irir scattering amplitude, where they 
represent the potentially largest two loop contribution. It is shown that their 
correction is reasonably small. 
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1 Introduction 



The tttt scattering reaction is of fundamental importance for understanding low en- 
ergy hadronic physics. Since the sixties the most fruitful way to look at it has been to 
consider the constraints given by the symmetry properties of the strong interactions 
Hamiltonian. The first predictions for the tttt scattering amplitude according to this 
method were worked out by Weinberg in 1966 ||]. In recent years all the theoretical 
work that goes under the name of Current Algebra has been incorporated in a more 
systematic and rigorous framework called Chiral Perturbation Theory (CHPT) B. 
Within this approach one assumes that the pions are the (pseudo) Goldstone bosons 
of QCD, and that the singularities generated by their exchanges dominate the Green 
functions at low energy After taking into account these singularities, CHPT works 
as a systematic expansion of their residues in powers of momenta and quark masses. 
Current Algebra corresponds to the leading term in this expansion. 

Gasser and Leutwyler have fully worked out the theory up to the one loop level, 
and in particular have determined all the new constants occurring in the effective 
Lagrangian at this level by comparison with experimental data || |J. After this 
determination one may calculate corrections to any Current Algebra prediction. For 
mi scattering this was done by the same authors || who found rather large corrections 
to Weinberg's predictions, even near threshold. It is interesting to try to understand 
why this happens. In fact, if we stick to threshold, we have only one expansion 
parameter, the average of u and d quark masses, which is known to be very small in 
comparison with the typical QCD scale. After expressing this parameter through the 
mass and decay constant of the pion it turns out to be (M^/AtiF^) 2 ~ 0.01. With 
such a small expansion parameter one would expect the corrections to the leading 
term to be rather small. However, the S-wave, 1 = scattering length, for example, 
gets at next order a 28% correction to the Weinberg prediction a[] = 0.16. While at 
first sight this may look surprising, the reason for having such large corrections is well 
known ||. Beyond leading order the unitarity property requires the appearance of 
nonanalytic functions such as M% logM^, which, if M 2 is small, may be significantly 
larger. In fact, for a[], 90% of the correction is given by the term proportional to 
Ml log M^/fi 2 , at a scale /i = lGeV. 

A natural question then arises, whether higher orders which contribute terms like 
(M^logM^)™ may still give large corrections. This is especially interesting because 
the comparison with experimental data is still not completely satisfactory, as can be 
seen in Table |l| For a{] in particular, large corrections would still be needed in order 
to obtain agreement with the present experimental central value, while Gasser and 
Leutwyler assigned a very small uncertainty (of the order of 5%, if one considers only 
the uncertainty in the value of the low energy constants) to their one loop calculation. 

On the other hand Stern and collaborators 0] have shown that by allowing the 
quark condensate to take an unexpectedly small value, and modifying accordingly 
the chiral expansion || , one could find agreement with the measurement of Rosselet 
et al. of a{] already at one loop (see Table 1 and Figure 2 in Ref. J7|). This makes 
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this quantity even more interesting, since the possibility for the quark condensate to 
be numerically near to F% (and not much larger) has not yet been clearly excluded, 
and this would be one of the few places where one could find experimental evidence 
against or in favour of it. 

Our knowledge about this quantity will most probably improve both on the ex- 
perimental and theoretical side in the near future. A full two loop calculation is now 
in progress |TIJ, and a new, high statistics measurement of decays which will be 



made at DA$NE, should sizeably reduce the present experimental error |TT . 

While waiting for the full two loop calculation, which, incidentally, is rather long 
and tedious, there is the possibility to evaluate at a rather low cost the potentially 
more dangerous part of the order 0(p 6 ), i.e. the double chiral logs. The calculation 
can be done in a relatively simple way by using equations that follow from general 
properties of the renormalization procedure. The content of these equations is exactly 
the same as that of the Renormalization Group Equations (RGE) of a renormaliz- 
able field theory, i.e. that (considering for example dimensional regularization) the 
residues of the double pole of two loop graphs are proportional to the residues of the 
single pole of one loop graphs. 

The use of these equations with a nonrenormalizable Lagrangian is not new. 



Weinberg derived them for the chiral effective Lagrangian in 1979 [12], while several 



authors have used them (under the name of pole equations) in the framework of the 
two dimensional a-model up to the four or five loop level p 



2 Renormalization Group Equations in CHPT 

As far as we know, the only application of RGE with an effective Lagrangian since 



Weinberg's original derivation |p!2| , has been Ref. ||14|| , where these equations have 
been used as a check on a full two loop calculation. 

Before deriving the equations we shortly review the basic ideas and notation of 
CHPT. Consider the QCD Lagrangian with two flavours in the isospin symmetry 
limit m u = = m, with external fields coupled to quark bilinears: 

£ = £qcd + <?7 m 0m + «M7s)g - q(s - i^p)q , (2.1) 

where £q CD is the QCD Lagrangian with massless quarks - the masses can be ab- 
sorbed in the external field s. 

Under the assumption that the pions, because of their small mass, dominate 
the low energy phenomena described by this Lagrangian, one can construct a good 
representation of the generating functional of QCD at low energy by means of an 
effective Lagrangian: 

J [dq][dA}e iC = e lZ[mp] = J [dU}e iC ^^ , (2.2) 

where [dA] stands for the integration over the gluon fields. The effective Lagrangian 
must be constructed with the only requirements of being the most general symmetric 
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Lagrangian containing the pions (through the matrix U) and the external fields. It 
can be expanded in a series of terms with increasing number of derivatives and powers 
of quark masses: 

C eff = C 2 + C 4 + C 6 + ... . (2.3) 

Each term L 2n contains a number of monomials with a coefficient which is not fixed 
by symmetry constraints. At present the best way to determine their value is by 
comparison with experimental data. C 2 contains only two arbitrary parameters, F 
and B, which are easily related to the decay constant and mass of the pion: 

F n = F(l + 0(m)), M 2 = M 2 (l + 0(rh)) , (2.4) 

where M 2 = 2Bm. At next order there are seven new free parameters. The relation 
between these constants called /j's, i = 1, . . . , 7, and measurable quantities has been 
worked out by Gasser and Leutwyler in Ref. ||. 

This expansion of the Lagrangian generates a corresponding series for the gener- 
ating functional: 

Z[v, a, s,p] = Z 2 + Z A + Z 6 + . . . . (2.5) 



As Weinberg has shown |12] , the contribution of loops to the generating functional is 
suppressed with respect to tree diagrams: Z 2 is given only by tree diagrams from £ 2 ; 
■Z4 receives contributions from tree diagrams from £4 and one loop diagrams with C 2 ; 
Zq contains tree diagrams from £ 6 , one loop diagrams with one vertex from £ 4 , and 
two loop diagrams containing only C 2 vertices; etc. For more details about CHPT 
we refer the reader to the original article by Gasser and Leutwyler || and to recent 
reviews [|]. 

Here we want to apply the RGE to the tttt scattering amplitude obtained from 
Zq. We derive the equations in a way which is much similar to the one followed in 



Ref. fll4[ . We use dimensional regularization and the minimal subtraction scheme. 
Before renormalization {i.e. before including the contribution of tree diagrams from 
£ 6 ), the 0{jP) 7T7T scattering amplitude looks like this: 



M 2t G{s, t, u) + M e J2 l iHi(s, t, 



(2.6) 



where e = d — 4. In this expression we have simply made evident the distinction 
between real two loop graphs, all incorporated into the function G(s,t,u), and the 
one loop graphs with one vertex from the £4 Lagrangian, that are included into the 
functions Hi(s,t,u). The dimension of the product F 6 A^ is (mass) 6+2e , then, to 
carry out the renormalization program we write 

= ft- {Laurent series of /i" 2e F 6 A (6) in e} . (2.7) 

The curly brackets contain a function with dimension (mass) 6 , that will be provided 
by powers of external momenta and masses. The poles in e appearing in the Laurent 
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series will be cancelled by the contribution of tree graphs from Cq. The ingredients 
we need to analyze the structure of the two loop amplitude are the following: 

nt + \ 92(s,t,u) gi(s,t,u) . 

G(s,t,u) = 5 1 \-G(s,t,u) + U(e) , 

e z e 

Hl (s,t,u) = Mf^f) + H i ( 3 ,t,u) + +0(e) , 

k = ^^ + lT s + 0(e^j (2.8) 

Note that G and the ifj's do not depend at all on /z, while the constants Zj's are by 
definition /i-independent; i.e. we have the equations: 

= Si , ^Tx = -l? S • (2-9) 
From general arguments on the structure of the divergences of loop graphs (see e.g. 



Ref. [fl5|), we know that the residues g 2 and h^s are polynomials in external momenta 
and masses. On the other hand, g\ is in general a nonlocal function of the same 
variables, due to the presence of divergent subgraphs in the two loop graphs. Despite 
the fact that CHPT is a nonrenormalizable theory, the renormalization program 
can be carried out in complete analogy with the case of a renormalizable theory. 
In particular, all the poles in e must have residues which are polynomials in the 
external momenta and masses, in such a way that one can remove them by adding 
counterterms to the original Lagrangian. If one writes down explicitly the Laurent 
series, one immediately realizes that for the double pole this is already obvious. In 
order to satisfy this property in the case of the single pole, the nonlocal functions 
appearing inside its residue must cancel. In other words, the following two equations 
must be satisfied: 

g 2 (s,t,u) = --^2Sihi(s,t,u) , 
1 i 

gi(s,t,u) = — 5jHj(s, t, u) + polynomial in s, t, u and M 2 . (2.10) 

i 

As a consequence of the first equation the residue of the double pole, and the coeffi- 
cient of the double log in the finite part, turn out to be proportional to the polynomial 

Ei 5ihi(s,t,u). 

With a relatively simple calculation we are then able to get interesting information 
about the two loops. To have the complete 0(p 6 ) result one still has to calculate 
part of the coefficient of the single log, and the finite contribution. 



3 Double logs in the tttt scattering amplitude 

We have calculated the polynomials hi partly using Feynman diagrams and partly 
with the Heat Kernel technique. In some cases we have performed the calculation 
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with both methods in order to have a cross check. More details about the calculation 
can be found in Ref. [16|. Here we simply give the results: 



hi(s, t, u) 



h 2 (s,t,u) 



h 3 (s,t,u) 



16tt 2 



16tt 2 



16tt 2 



_^l s 3 _ 5 / _ 2 200 M2 2 _ 4 2 _ 2 
2 6 v ; 3 3 v ; 

_164M 4 S + — M 6 
3 



107 



25 , 
12 v 



12 

296 „ , 208 , - 

M 4 s + M 

3 3 



-4M 4 s - 2M e 



44MV - -M 2 (t-u) 2 
3 V ; 



(3.11) 



All the other h^s are zero. Before using these expressions to get numerical results, we 
need to shuffle part of the two loop amplitude into the lower order amplitudes 
in order to express everything in terms of the physical pion mass and the physical 
decay constant. Since we are calculating only the divergent part of one loop graphs 
containing one U in we need to calculate the same contribution to M 2 and F w : 



Ml 



M 



' M 2 



2/ 3 + 



1 A(0) 



1 + 



+ 



M 4 



M 2 
M 4 



2L - 



2i M 2 
14/i - 8/ 2 - 
2A(0)" 



3/) 1A(0) + 



+ 0(M 6 



F 4 - 



i M 2 

(14^ + 8/2-4/3-3/4)- 



1 A(0) 

i M 2 



+ 0(M 6 ) ,(3.12) 



where A(z) is the Feynman propagator of a scalar field with mass M: A(0) = 
2iM7(167r 2 e) + 0(1). 

In order to evaluate the size of these corrections to the nir scattering amplitude, 
we calculate their contribution to threshold parameters, following the conventions of 
Ref. IH 1 : 



7M 2 
32ttF 2 

1 

W 2 



1 + 



1 + 



+ 



Mi 

7T 

IT 

Mi 
~F 4 



58 

T 



96 

y 



-ko - 5fe 



11 



86, 73, 5, 29, 

Ki Ko K-i Ka + 

3 2 2 3 



+ 0{Mt) 



+ 0{Mt) 



lr £\ie ellipses between parentheses stand for the one loop contributions that we do not display 
here in order not to put too many formulae in this note. We refer the reader again to Ref. || (note 
that the corrections there arc usually expressed in terms of M and F and not with the physical 
mass and decay constant, as we do here). 
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M 2 { M 4 
— + ..) + — 

8ttF2 I 1 ^ 

7T ^ 7T 

1 f M 4 



2/ci + 6/c 2 + /c 3 + -/c 4 + . . . 

20, 3, 5, 7, 
— 3~ + 2 _ 2 + 3 + 



119 5 

■12fci — k 2 - -h - 3k 4 + 

o 2 



18ttF6 



-23fci - — k 2 + ■■■ 



+ 0{Mt) 




«2 



1 

Ft 



1 

F 4 ' 



•) + 



Ml 



lOvrFf 



IOttF 6 



11. 8, 5, 

— k\ ko kd + . . . 

18 3 9 

16, 3, 2, 

— k\ H — fco + . . . 

9 4 9 



+ 0{Mt) 



+ 0(Mi) 



(3.13) 



where 



167T : 



, M l 
/i 2 



(3.14) 



Note that with this definition of the fcj's we are also taking into account in Eqs. 
( p. 13 ) the part of the single log which is proportional to the renormalized constants. 
From our previous considerations one can see that this contribution is unambiguously 
determined by the polynomials Note also that the renormalized constants we 
are using here are defined in the renormalization scheme adopted by Gasser and 
Leutwyler in Ref . , and not in the MS scheme that we used in the previous section 
to simplify the notation. 

Before analyzing the part of the two loop contribution evaluated here it is in- 
teresting to investigate in some more details the structure of the 0(p 4 ) correction. 
In the third column of Table [l[ we show how much of that correction is due to the 
chiral log evaluated at a scale \i = lGeV. It is clear that for the S'-wave parameters 
the chiral log is responsible for almost all the correction. On the other hand, for the 
P-wave, the chiral log arises only through F n renormalization - had one used the 
pion decay constant in the chiral limit F, there would not have been any log (and 
in fact for b\, which does not contain any tree level contribution, so that the use of 
F or F n is equivalent, there is no chiral log at this order). Moreover, in this case 
we expect the p to give an important contribution, and it is well known that the 
resonances manifest themselves through the Z['s. Hence the fact that the chiral logs 
are not the main part of the one loop correction does not come out as a surprise. 
Finally we have the two D-w&ve scattering lengths, for which the order 0(p 6 ) is the 
next to leading correction. In this case we feel that it is more hazardous to estimate 
it only through the double chiral logs; also because, a posteriori, the part of the single 
chiral logs that we get turns out to be as important as the double logs. Moreover 
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the 1 = 2 scattering length at leading order is unnaturally small because of a strong 
cancellation between the chiral log and the finite part. 

From all these observations we conclude that for P and D waves the assumption 
that the double chiral logs could be a good estimate of the full two loop correction 
is poorly justified. 

The numerical results are displayed in Table [1], together with the one loop pre- 
dictions and experimental data available. Since the numbers we produce are in- 
dependent we have calculated them for two reasonable values of \i. Moreover, in 
Table [l] we show both the double logs and the complete fcj's contributions. 

We observe that: 

1. For the S- wave parameters this 0(p 6 ) correction is, in the extreme case, of the 
order of 10%. To judge whether the size of this correction is reasonable or not, 
we may compare the coefficients of the single and double logs in the 0(p A ) and 
0(p 6 ) corrections, respectively: 



where L = (M 2 /16n 2 F 2 ) log M^//i 2 . In first place one can see that the ab- 
solute value of the L 2 coefficient is numerically quite similar to the square of 
the L coefficient, in all cases but 6q. [Incidentally, had respected this rule 
of thumb for the absolute value, with the sign it has, the distance from the 
present experimental central value would have increased sizeably.] The second 
interesting point is that they show a very regular sign pattern. In order to have 
some more insight into it, it is useful to disentangle two effects of comparable 
importance: the correction to the amplitude A(s,t,u) expressed in terms of F, 
and the effect of the renormalization of F into F n . This is done in Table |3], 
where one can see that whilst the two effects add up for the 1 = parameters, 
they partially cancel in the case of the 1 = 2 parameters. In particular, in the 
latter case one can trace back the minus sign of the coefficient of L 2 to the 
interference of the two one loop corrections. 

2. As far as the P- and D-wave parameters are concerned, we have displayed here 
the double chiral logs contribution, essentially for the sake of completeness. As 
we stressed already, we see no special reasons why the numbers we find here 
should be representative of the full two loop correction. 



o 
o 







(3.15) 
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As we mentioned in the introduction, it is hoped that there will be an im- 
provement in the experimental knowledge of a[] in the near future, arising from 
a new measurement of K e 4 decays. On the theoretical side the full two loop 



calculation, which is now in progress l lQfl , will improve the one loop result of 



Gasser and Leutwyler |3], [5|]. Our results show that the contribution from the 
two loops is expected to be reasonably small, and confirm that it looks very 
unlikely that CHPT could at any order reproduce the present experimental 
central value. 
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Table 1: Numerical values of the threshold parameters according to Eqs. ( |3.13| ), for 
\x = 0.5, lGeV. In the column containing the one loop results, the numbers between 
parentheses show how much of the one loop contribution is due to the chiral log 
evaluated at the scale \i = lGeV. In the calculation we have used M n = 139.6MeV 
and F n = 93.2MeV. To calculate the fcj's we have used the Z['s given in Table 
Scattering lengths aj, and effective ranges b\ are given in the appropriate powers of 
M w +. The value of /i is expressed in GeV. 





1 loop (log) 

M = l 


1 loop +log^ Ml 
\i = 0.5 [A — 1 


1 loop+Zcj's 
fi = 0.5 [A — 1 


experiment 
17 


a 

-10 a 2 
-10 bl 
10a\ 


0.201 (88%) 
0.248 (125%) 
0.418 (133%) 
0.726 (100%) 
0.370 (46%) 


0.205 0.211 
0.260 0.277 
0.415 0.411 
0.721 0.713 
0.379 0.390 


0.211 0.213 
0.275 0.279 
0.409 0.407 
0.701 0.691 
0.391 0.395 


0.26 ±0.05 
0.25 ±0.03 
0.28 ±0.12 
0.82 ±0.08 
0.38 ±0.02 


10 2 b\ 
10 2 a° 2 
10 3 a2 


0.485 (0%) 
0.181 (122%) 
0.205 (432%) 


0.619 0.805 
0.206 0.241 
0.145 0.061 


0.749 0.780 
0.273 0.302 
0.267 0.335 


0.17 ±0.03 
0.13 ±0.3 
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Table 2: Numerical values of the coupling constants l^s and the k^s used in Table 
[I]. The values of /g and l\ are taken from Ref. 0], while for the other two we have 
used the more recent ones obtained in Ref. |18|| . The value of /i is expressed in GeV. 



i 






lO 3 ^ 




/i = 0.5 /i 


= 1 


/i = 0.5 fi = 1 


1 


-4.5 


-5.9 


0.20 0.39 


2 


7.5 


4.6 


-0.65 -0.87 


3 


-0.5 


1.6 


0.17 0.15 


4 


11.1 


2.3 


-1.24 -1.47 



Table 3: Breakdown of the coefficient of L and L 2 into contributions arising from two 
different effects: the corrections to the amplitude A(s,t,u), and the renormalization 
of F into F n . For the L 2 coefficient we have also separated out the contribution due 
to the interference between the two one loop effects. [The small effect due to mass 
renormalization is included into the columns tagged by A(s,t,u).] 





L-coe 


fficient 


L 2 -coefficient 




A(s, t, u) 


F n 


total 


A(s, t, u) 


interference 




total 


"o 


5 
~2 


-2 


9 
~2 


145 
21 


10 


7 
2 


857 
42 




20 


-2 


26 


785 


80 


7 


1871 








36 


y 


2 


36 


a 2 


7 
2 


-2 


3 
2 


16 

y 


-14 


7 
2 


31 

~y 


b 2 


16 


-2 


10 


476 


64 


7 


83 




~3~ 




y 


36 


~y 


2 


~18 
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